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Abstract 

The generation of harmonic radiation through a non-hnear mechanism driven by 
bunching at fundamental frequency is an important option in the operation of high 
gain Free-Electron Lasers (FELs) . The use of harmonic generation at a large scale 
facility may result in achieving shorter radiation wavelengths for the same electron 
beam energy. This paper describes a theory of second harmonic generation in planar 
undulators with particular attention to X-Ray FELs (XFELs). Our study is based 
on an exact analytical solution of Maxwell equations, derived with the help of the 
Green's function method. On the contrary, up-to-date theoretical understanding 
of the second harmonic generation is only limited to some estimation of the total 
radiation power based on the source part of the wave equation. Moreover, we find 
that such part of the wave equation is presented with several incorrect manipulations 
among which is the omission of an important contribution. Our work yields correct 
parametric dependencies and specific predictions of additional properties such as 
polarization, angular distribution of the radiation intensity and total power. The 
most surprising prediction is the presence of a vertically polarized part of the second 
harmonic radiation, whereas up-to-date understanding assumes that the field is 
horizontally polarized. Altogether, this paper presents the first correct theory of 
second harmonic generation for high gain FELs. 
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1 Introduction 



In a Free-Electron Laser the electromagnetic field at the fundamental har- 
monic interacts with the electron beam. As a result, the beam is bunched in 
non-linear (sinusoidal) ponderomotive potential. When the bunching is strong 
enough, the beam current exhibits non-negligible Fourier components at har- 
monics of the fundamental as well. In the SASE case only the transverse 
ground mode of the fundamental harmonic survives, due to the transverse 
mode selection mechanism P in the high gain regime, and is responsible for 
the bunching mechanism. As a result, the nonlinear Fourier components ra- 
diate coherently and the phenomenon is referred to as (nonlinear) harmonic 
generation of coherent radiation. 

The process of harmonic generation of coherent radiation can be considered 
as a purely electrodynamical one. In fact, the harmonics of the electron beam 
density are driven by the electromagnetic field at the fundamental frequency, 
but the bunching contribution due to the interaction of the electron beam with 
the radiation at higher harmonics can be neglected. This leads to important 
simplifications. In fact, in order to perform numerical analysis of the char- 
acteristics of higher harmonics radiation, one has to solve the self-consistent 
problem for the fundamental harmonic only. Subsequently, the solution to this 
problem, that must be obtained with the help of a self-consistent code, can be 
used to calculate the harmonic contents of the beam current. These contents 
enter as known sources in our electrodynamical process: solving Maxwell equa- 
tions accounting for these sources gives the desired characteristics of higher 
harmonics radiation. As a result, simulation codes dealing with harmonic gen- 
eration are not first principle codes: in fact, they simply compute the solution 
of Maxwell equations obtaining the proper sources by means of first-principle 
codes. 

Non-linear generation of the second harmonic radiation, in particular, is im- 
portant for extending the attainable frequency range of an XFEL facility: in 
fact, the high peak-brilliance increase of XFELs with respect to third gen- 
eration light sources (up to eight orders of magnitude) makes the second 
harmonic contents of the XFEL radiation very attractive from a practical 
viewpoint. Moreover, it is also important in connection with experiments that 
make use of the fundamental harmonic only. In fact one must be able to esti- 
mate correctly the higher harmonics effects to distinguish between nonlinear 
phenomena induced by the fundamental and linear phenomena due to the 
second harmonic. 

The subject has been a matter of theoretical studies in a high-gain SASE FEL 
both for odd and even harmonics 0,0,0]; where the electrodynamical prob- 
lem is dealt with. The practical interest of these studies is well underlined by 
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the fact that they were followed by both numerical analysis [5| and experi- 
ments, that have been carried out in the infra-red and in the visible range of 
the electromagnetic spectrum Q . 



Experimental results are compared with numerical analysis and numerical 
analysis rely on analytical studies: this fact stresses the importance of a cor- 
rect theoretical understanding of the subject. Remarkably, such understanding 
does not require the introduction of radically new physical mechanisms. The 
key ideas involved are not much different from those regarding the second har- 
monic generation from a single particle, treated a long time ago and presented 
in Synchrotron Radiation textbooks (e.g. llOl)- A complexity though, is 
constituted by the presence of many electrons involved in the radiative pro- 
cess, each with a given offset and phase, radiating coherently, as a whole or 
in part, due to the longitudinal modulation of the beam current at the second 
harmonic. 

For a given frequency component the electromagnetic wave equation dic- 
tates both a characteristic longitudinal length (that is the radiation formation 
length) and a characteristic transverse length. As we will see, when the beam 
transverse size is smaller than the characteristic transverse length the entire 
electron beam behaves like a single electron and the harmonics of the beam 
current are simply interpretable in terms of the harmonic contents of a single 
particle current: in this case, all the particles act coherently and the radiated 
intensity scales with the square number of the electrons in the beam. When 
the transverse size of the beam increases and becomes much larger than the 
characteristic transverse length less electrons contribute collectively to the 
field and if the beam current remains constant the total radiated power is 
decreased. 

The characteristic transverse length is specified in a natural way after a di- 
mensional analysis of the problem. The first treatment 0] of non-linear gener- 
ation of even harmonics does not account for the presence of such parameter, 
followed in this by others 0, S] . We find that these works include arbitrary ma- 
nipulations of the source terms in the paraxial wave equation. Among these, 
an important part of the source terms is systematically dropped. Moreover, 
estimations of the second harmonic power are based on the electromagnetic 
sources (after manipulation) while exact calculations should be based on a 
solution of Maxwell equations. Altogether, we find that these works predict 
an incorrect dependence of the second harmonic field on the problem parame- 
ters. Results of 0] are extended in to the case of an electron beam moving 
off-axis through the undulator. One of the conclusions in [8] is that the second 
harmonic power increases when an angle between the beam and the undula- 
tor axis is present. We find that the power of the second harmonic radiation 
should never increase when such angle is present: in particular, as we will see, 
it is independent of it in optimal situations when the microbunching wavefront 
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is matched with the beam propagation. 



In this paper, that was inspired by a method developed to deal with 
Synchrotron Radiation from complex setups, we present a theory of second 
harmonic generation in high-gain FELs. First we give, in Section 2, an exact 
analytical solution of the wave equation for the second harmonic generation 
problem. The procedure employed to derive such a solution shows the advan- 
tages of a Green's function method. In Section 3, our result is used to calculate, 
in a particular case, specific properties of the second harmonic radiation such 
as polarization, directivity diagram and total power including proper para- 
metric dependencies. The most surprising prediction of our theory is that the 
electric field is not only horizontally polarized, as it is usually assumed, but 
exhibits, though remaining linearly polarized, a vertically polarized compo- 
nent too. Following the presentation of our theory, in Section 4 we comment 
on the differences between our approach and the present understanding of 
the second harmonic generation mechanism. Finally, in Section 5, we come to 
conclusions. 



2 Complete analysis of Second Harmonic Generation mechanism 



As has been said in the Introduction, the process of (second) harmonic gen- 
eration of coherent radiation is a purely electrodynamical one. First, proper 
initial conditions are given as input to an FEL self-consistent code, which 
calculates the electron beam bunching from the interaction of the beam with 
the first harmonic radiation. Then, the results from the self-consistent code 
are used as electromagnetic sources to solve the problem of second harmonic 
generation. For simplicity, in the following we will consider a beam modulated 
at a single frequency uo as the source. One may always write the longitudinal 
current density jz along the undulator as a sum of an unperturbed part inde- 
pendent of the modulation and of the time, jo, and a term responsible for the 
beam modulation, j^, at frequency u (perturbation): 

3z{z,f±,t) = joiz,f_i) + jz{z,f±,t) . (1) 



We assume that we can write the unperturbed part jo as if all the particles 
where moving coherently, that is 

jo{z, rl) = jo{r± - rf\z)) , (2) 



where ry (^) describes the coherent motion. This assumption is always veri- 
fied, for instance, in the case of a single particle, when jo is simply a (5-Dirac 
function, or in the case of a monochromatic beam. If some energy spread is 
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present, in order for Eq. (2) to be valid we should assTime that the transverse 
size of the electron beam is not smaller than the typical wiggling motion of 
the electrons. In this case, the validity of Eq. (2) has an accuracy given by the 
relative deviation of the particles energy form the average value, 57/7. Since 
for the PEL process ^7/7 is, at most, of the order of the efficiency parameter, 
we have ^7/7 -C 1 and Eq. (2) is valid with the same accuracy of EEL theory. 
However it should be noted here that the average energy of the beam is to 
be considered, in general, a function of the coordinate 2;, 7 = 7(2;): it has to 
be given as a result of start-to-end simulations and considered as an input for 
our electromagnetic problem. 

The perturbation can then be written as 



X a2 {z, r_L — rf^{z) ) exp 



f dz' 



+ C.C. 



(3) 



The function 02 is to be considered a result from the EEL self-consistent code, 
and its dependence on z describes the evolution of the modulation through 
the beamline and accounts for emittance and energy spread effects. It should 
be noted that the values of 0,2 are not necessarily real: in fact there can be a 
z-dependent phase shift with respect to the phase uj /q dz' /vz{z') — cut. 

In order to correctly calculate the phase uj Jq dz'/vz{z') — cot in Eq. (3) one has 
to account for the dependence of the longitudinal velocity associated with the 
coherent motion, Vz on the position z. The function Vz{z) can be recovered 
from the knowledge of r^f^{z) and of the average energy of the beam 7 = ^{z). 

If the beam is deflected of angles rjx and rjy in the horizontal and vertical di- 
rection with respect to the z axis, the velocity of the coherent motion depends 
also on the deflection angles. Renaming position and velocity of the coherent 
motion with no deflection with the subscript "(nd)" one obtains: 



1 - J 

V±{Z: V) = V±(nd) (z) + V^i^d) {z)f] , 



(4) 



and 



rf^ {z,fj)^ fi'^ (z) +fjz . (5) 

Also 0,2 will depend on ff. The exact dependence is flxed by the way the beam 
is prepared and should be regarded as a condition for the orientation of the 
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microbunching wavefront. In all generality we can write: 



(6) 



In the limit for 7^ ^ 1, the total current density can be written as 



^..(rx-rt'(.,.-))x(l + 



^2 (z,fx_ - r-l\z,f])^ 



exp 



dz' 



— iujt 



u 

One can express the charge density as 



+ C.C. 



(7) 



Jz ]z 

p = — ^ — , 

V,, C 



(8) 



as we will be working in the paraxial approximation. 



Eq. (7) and Eq. (8) give us the expressions to be used as sources for Maxwell 
equation. Looking for solutions for Ex in the form 



E± = Ei_ exp [iuj{z/c - t)] + C.C. 



(9) 



and applying the paraxial approximation, one may write the Maxwell equation 
describing Ej_ as 
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V±' + — E^ = — exp 

, c oz I c 







'iuj ^ -1 ' 















joa2 , (10) 



were we have put 



dz' 



Vzi.z'.ff) 



:iii 



With the aid of the appropriate Green's function an exact solution of Eq. 
(10) can be found without any extra assumption about the parameters of the 
problem. 



~ 11 
E±{zo,r^o) = — / dz' / dr'^ 

C J Zo — Z' J 



-rV±{z',r])-V^ 



7 



xjo (r\ - ff(z', fj}) ~a, (z\ r\ - 7^l\z' , ff)) 



exp 



r_Lo - r 



I. |2 



2c{zo - z') 



^s{z\ff) -UJ- 



(12) 



where V'j_ represents the gradient operator with respect to the source point, 
while {zo,'r±o) indicates the observation point. Integration by parts of the 
gradient terms leads to 



E^^ — 2 / 

C"^ J Zo — Z' 

—oo 



dr 



) (v±{z',fi) f±o-r'j 



z„ - Z' 



xjo {r'± - rf\z', fj)^ 0,2 [z', r'± - rl\z', 77)) exp i^riz', r'±, fj) 



(13) 



where the total phase $t is given by 

^ . |2 



C 



+ 0; 



\r±o - r 
2c{zo - z') 



(14) 



We will now make use of a new integration variable I — r'± — r"^ [z' , ff) so that 



El = 



and 



'dz'- ' 



Zo - Z' 



dl 



Vl_{z' , f)) f±o - f^±\z', ff) - I 



Zn - Z' 



xjo (l) 0.2 {z\ I) exp [i$T(2;', I, ff) 



f^±o - rf^jz^ff) - l\ 
2c{zo - z') 



(15) 



(16) 



We will consider the case of a planar undulator and we will be interested in the 
total power of the second harmonic emission and in the directivity diagram 
of the radiation in the far zone. Accounting for the beam deflection angles 77-j: 
and r\y we model the electron transverse motion as: 



vi_{z',ff) 



sm {kyjZ ) + ri^v^ 

7 



^ + [VyVz] y , 



(17) 



and 
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K 

iky, 



(cos {kyjz') - 1) + r]xz' + 



x+[vyz' + ly]y . (18) 



Here K is the deflection parameter and ky, = 27r/A^, being the undulator 
period. Moreover, one has 



4^2 



^ Vx+Vy \ ^ _ KV: 



kwl 



^^^^jin {2k^z') + -^ cos ik^z') 



(19) 



We will now introduce the far zone approximation. Substitution of Eq. (19), 
Eq. (18) and Eq. (17) in Eq. (13) yields the field contribution calculated along 
the undulator: 



E± = J dl j dz'jo (f) d2 [z', f) exp 

-Lu; /2 



Y sin {ky,z') + {9^ - r]x)j x + {9y - %) y 



(20) 



where 



+UJ 



^^+^+i'{io^-v.r+{0y-vyr) 

sin {2kwz') — cos {kyjz') 1 

ikyjC J 

+ 9yly) + {91 + 91) 



^I'^ky.C 

K 



kuijc c 



2c 



(21) 



Here 9^ and 9y indicate the observation angles Xq/zq and yo/zo- Moreover, the 
integration is performed in from — to Ly,/2 in dz', = N^X^ being the 
undulator length. In fact, working under the resonance approximation in the 
limit for A^^„ 3> 1 allows us to neglect contributions outside the undulator [ll| . 



We will make use of the well-known expansion (see |l2l |) 

oo 

exp[iasm (ip)] = ^ Jp{a) exp [ipip] , 



f22) 



p=— oo 



where Jp indicates the Bessel function of the first kind of order n. 
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We will be interested in frequencies around the second harmonic: 

i^2o = ^kyjC^l , (23) 



where 



1 + 7^2/2 



(24) 



Indicating with E_i^2 the second harmonic contribution calculated at frequen- 
cies around ui2o one obtains 



. oo oo 



dz'jo (f) a2 f) exp[i^o] 



-Lyj/2 

^ J2 J2 Jmiu)Jniv)exp 



—oo — oo 

oo oo 



m=— oo n=— oo 



zrni 



X< 



+{^x - Vx) exp {iR^ky^z'} 



x + 



(^2/ - Vy) exp {iRojkyjz'} 



where 



UJ 



Ruj = n — 2m , 



y\ , (25) 



(26) 



with 



-1 



(27) 



Moreover 



u 



[1 - irV(47')] 



^1 4 {l + f + 72 [{6, - r^^f + [Oy - r^yf] } ' 

^20 2K- [l-A'V(lvO](^.r -'/.-) 



^ll + f +7^ (^. - ^x)' + (^. - ^.)^ 



(28) 
(29) 



and 



A;^7c 



{e. - vx) - -{Oxix + eyiy) + Yj.oi + e'y) 



(30) 
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In the limit for A^^ ^ 1 and if 02 does not vary much in z' over a period of the 
undulator A^^, the fast oscillations in the exponential function in the integrand 
of Eq. (25) tend to suppress the integral unless = 0, = —1 or R^ = 1, 
that is when at least one of the exponential function is simply unity. U u = 002 
this corresponds to n = 2 — 2m, n = 3 — 2m and n = 1 — 2m respectively. 
Neglecting all other terms and imposing uj = 002 + Auj2 we obtain 



E 



±2 ■ 



iUJ2o 
C^Zo 



00 



00 


Lw /2 




/ 


—00 




.K 






( 


"'2^ 





dz' jo ipj a2 (^z' , exp[i$o] exp 
[ i[3 — 2m]7r 1 



. Ato'2 



1/ ^ 

UJi 



-Jra{u)Jl-2m{v) CXp 



i[l — 2m]7r 



+ - Vx)Jm{u)J2-2m{v) CXp 



+ 



i^y - Vy)Jm{u)J2-2ra{v) CXp 



' i[2 - 2m]7r| 
■i[2 -2m]7r' 



X 



(31) 



For any value of K and Ox — f]x much smaller than 1/7^ , f is a small parameter 
and only the smallest indexes in the Bessel functions Jq{v) ~ f ^ in Eq. (25) give 
non negligible contribution. As a result, Eq. (25) can be drastically simplified. 
One can write (compare, for instance, with |l l| ) : 



00 00 



E 



±2 ■ 



l^2o 
C^Zo 



dz' exp [i$o] 



-00 —00 



j2 



xjo (^') 02 (^',/) exp[iCz']exp|i27^ {Ox - rj-xf + {Oy - rjyf k^^z'^ 

(32) 



where we have defined 



A 



2K' 



B = Ji 



2 + ^2 



2 + K^ 



2 + K^ 



-J 2 



2 + K^ 



+ Ji 



2 + K^ ' 



(33) 
(34) 



we have used the fact that 
Aa;2 



H % - Vx) + {Oy - Vy) + c , 



(35) 
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Fig. 1. Illustration of the behavior of the ratio Q{K) between the second harmonic 
field contribution due to the gradient of the density part of the source and the 
contribution due to the current part of the source. 

with 



c 



^2o 



(36) 



and, under the resonant approximation, we have 



C 



-{ej^ + Oyiy) + -^{el 



(37) 



The detuning parameter C should indeed be considered as a function of 2, 
C = C{z) which can be retrieved from the knowledge of 7 = 7(^). 

It is important to see that the terms in Ji in Eq. (33) and Eq. (34) are due 
to the presence of the gradient term in V_l(jo«2) in Eq. (10), which has been 
omitted in [i^ and later on in We find that, without the gradient, term 

one would recover results quantitatively incorrect for the x-polarization com- 
ponent. In Fig. 1 we plotted the ratio between the contribution of the radiation 
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field due to the gradient of tlie density part of tlie source and the contribution 
due to the current part of the source for the x-polarization component. This 
is a function Q{K) of the K parameter only and it can be written as 



Q 



E 



±2g_ 

'±2c 



2 + 
2^2 



2 + 



2 + K\ 



-J2 



2 + K\ 



(38) 



where the subscript "g" stands for "gradient and "c" stands for "current". 
As it can be seen from Fig. 1, the gradient term always contributes for more 
than one fourth of the total field, independently of the values of K. Also, if 
the gradient term is omitted, the entire contribution to the field polarized in 
the y direction would go overlooked. The inclusion of the gradient term in the 
source part of the wave equation should not be considered as a peculiarity of 
the second harmonic generation mechanism. In Synchrotron Radiation theory 
from bending magnets, for instance, the presence of such a source term is 
customary and it is responsible, as here, for part of the horizontally polarized 
field and for the entire vertically polarized field. Moreover, the gradient term 
is always associated with an integration by part, and therefore is always ac- 
companied with the gradient of the Green's function, which is responsible for 
a term proportional to the observation angle Ox,y 

Eq. (32) can be also written as: 



E±2 = -IT- exp 



C^Zo 



00 00 



^ {f^x^x ~l~ ^y^y) 



X j dlx j dly j (i2;'exp 
—00 

'x-r]xf + {ey-r^yf] z'\p^^\z'XC) , 



-00 —00 



X exp <^ I — 
^ ^ 2c 



where we have defined p as 

p(2)(/, C) = Jo if) a2 (/, T) exp [iCz'] Hl^z') , 



(39) 



(40) 



Hl^{z') being a function equal to unity over the interval [—Lw/2, Lm/2] and 
zero everywhere else. Its introduction simply amounts to a notational change. 
Namely it accounts for the fact that the integral in dz' is performed over the 
undulator length in Eq. (32), while it is performed from —00 to 00 in Eq. (39). 
It should be noted that, usually, computer codes do not present the functions 
0,2 and exp[iC2;'] separately as we did, but rather they combine them in a 
single product, usually known as the complex amplitude of the electron beam 
modulation with respect to the phase = 2kyjz' + {u/c)z' — cot. Regarding p^"^^ 
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as a given function allows one not to bother about a particular presentation 
of the beam modulation. 



Eq. (32) or, equivalently, Eq. (39) are our most general result, and are valid 

independently on the model chosen for the current density and the modulation. 
It is interesting to note here that, when one writes Eq. (32) in the form of Eq. 
(39), one obtains an expression which is formally similar to the spatial Fourier 
transform of p^'^\z',l,C) with respect to z' and /. There are two problems 
though: first, p^^^ is a function of ff, which appears in the conjugate variable 
to z' and, second, if 7 = 7(2;) one has UJ20 — <^2o(-2^)- 



3 Analysis of a simple model 



Let us treat a particular case. Namely, let us consider the case when we can 
consider 7(2;) — ^ — const, when C{z) — and 



(^^^ ; j ^ (Jj exp 



(41) 



with 020 = const and 






(42) 



lo and a being the bunch current and transverse size respectively 

This particular case corresponds to a modulation wavefront perpendicular to 
the beam direction of motion. In this case Eq. (39) can be written as 



E±2 = — ^ exp 



C^Zn 



[A{e,-Vx)x + B{9y-riy)y\ 



00 00 



J dlx J dly J dz'exp |-^-^ [i^x - Vx) k + {dy - Vy) 



X 



X exp < I 



2c 



{Ox - Vxf + {Oy - Vyf] z' \ jo (f) HlM' 



(43) 



and amounts, indeed to the spatial Fourier transform of jo (f^ Hl^{z'). We 
obtain straightforwardly: 



iIoa2oU!2oLw 

E±2 = ^ exp 



C^Zn 



i^zoiel + el) 



[A{9x-rjx)x + B{9y-rjy)^ 
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xsmc 



X exp 



Lw^2o r 
. 4c ^' 

2c2 




(44) 



If the beam is prepared in a different way so that, for instance, the modulation 
wavefront is not orthogonal to the direction of propagation of the beam, Eq. 
(39) retains its validity. However it should be noted that, in this case, Eq. 
(39) is not, in general, a Fourier transform. It is if 7(2;) = 7 = const and 



. The case a = 1 has just 



p^^^ includes a phase factor of the form exp iafj ■ I 
been treated. The case a = corresponds, instead, to a modulation wavefront 
orthogonal to the z axis and not to the direction of propagation. 



Going back to our particular case in Eq. (44), a subject of particular interest 
is the angular distribution of the radiation intensity along the x and y po- 
larization directions which will be denoted with l2{x,y)- Upon introduction of 
normalized quantities: 




one obtains 



h{x,y) (4 -flx,6y- fjy) = COUSt X (fj^^y - 



xsmc 



X 



{Ox - Vx) + {Oy - fly) 

eXp|-A^ (^0^ - fl^y + (Oy - f]y 



(47) 



Note that in the limit for <^ 1, Eq. (47) restitutes the directivity diagram 
for the second harmonic radiation from a single particle. In agreement with 
Synchrotron Radiation textbooks 0| none of the polarization components 
of l2{x,y) has azimuthal symmetry, contrarily with what happens for the first 
harmonic, where only the x polarization is present and is endowed with az- 
imuthal symmetry. 
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Fig. 2. Plot of the directivity diagram for the radiation intensity as a function of 
0x — Vx Ski 9y — fjy = for the horizontal polarization component, for different values 
of A^. 

As an example, the directivity diagram in Eq. (47) is plotted in Fig. 2 for 
different values of A?" as a function of 9^ — fj^ Sii 9y — fjy — Q for the horizontal 
polarization component. 

The next step is the calculation of the second harmonic power. The power for 
the X- and y-polarization components of the second harmonic radiation are 
given by 



W2{x,y) = -^ j dxo J dyo\E±^^y{zo,Xo,yo,t)\''' 



—oo — oo 
oo oo 

c 



-oo — oo 



where (...) denotes averaging over a cycle of oscillation of the carrier wave. 

We will still consider the model specified by Eq. (41) and Eq. (42) with C = 0. 
It is convenient to present the expressions for and W2y in a dimension- 
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Fig. 3. Illustration of the behavior of F2{N). 

less form. After appropriate normalization they both are a function of one 
dimensionless parameter only: 



47V2 



(49) 



Here W2x = ^2x1^03^ W2y = W2y/W^y'' are the normalized powers, while 
the normalization constants W^^ and Wj^y^ are given by 



i3 



2 27rc 



(50) 



For practical purposes it is convenient to express Eq. (50) in the form: 



"■2o 

27r 



(51) 
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where Wb = nieC^'ylo/e is the total power of the electron beam and I a = 
rrieC^/e ~ 17 kA is the Alfven current. 

The function F2{N) is plotted in Fig. 3. The logarithmic divergence in F2{N) in 
the limit for A^ ^ 1 imposes a limit on the meaningful values of A^. On the one 
hand, the characteristic angle 6'max associated with the intensity distribution is 
given by 0^^^^ ~ 1/A^. On the other hand, the expansion of the Bessel function 
in Eq. (32) is valid only as 6"^ < A^^. As a result we find that Eq. (49) is valid 
only up to values of A^ such that A^ > N~^. However, in the case A^ < A''"-'^ we 
deal with a situation when the dimensionless problem parameter A^ is smaller 
than the accuracy of the resonance approximation ~ In this situation 

our electrodynamic description does not distinguish anymore between a beam 
with finite transverse size and a point-like particle and, for estimations, we 
should make the substitution In (A^) — > In (A^,^ 



' w ) 



We will now compare our results for the second harmonic with already known 
results for the first. The case treated in 3| corresponds to a modulation 
wavefront orthogonal to the direction of propagation, exactly as specified here 
for the second harmonic (i.e. perfect resonance with Eq. (41) and Eq. (42) 
valid) and allows direct comparison of results. The outcomes of have been 
presented, similarly to what has been done here for the second harmonic, in 
dimensionless form. After appropriate normalization, one finds: 



Wc 



(1) 



Fi{N) = - 

77 



arctan + - In 



(52) 



where the normalization factor Wjp is given by 



2 + ^2 



(53) 



Ajj being given by 



A 



jj 



Jo 



4: + 2K^ 



Jl 



4: + 2K^ 



(54) 



Here aio is the analogous of 020 for the first harmonic. For notational reasons, 
aio is one half of the original modulation level in Eq. (27) of jjl^]. It should 
also be noted that all A^ in Eq. (52) are multiplied by a factor 1/2 with respect 
to what is reported in [13]. This is because we are referring all results to the 
Fresnel number for the second harmonic. 

The function Fi{N) is plotted in Fig. 4. We can compare more quantitatively 
the normalized power for the second and for the first harmonic: 
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Fig. 4. Illustration of the behavior of Fi{N). 



F2iN) ^ vr 



1^ (i + im) 



arctan(^) + f ln(^)_ 



(55) 



while, from a practical viewpoint, the comparison between the real powers is 
equal to 



1 2 + K^ al A' + B^F2{N) 

F,{N) ~ (27r)37V, al A^jj F,{N) 



(56) 



It is interesting to calculate Eq. (55) in the limit ^ 1. We have 



F2(N) TT 



(57) 



In Fig. 5 we plot the behavior of F2/F1 as a function of N and its asymptotic, 
n/im), for A^> 1. 

Finally, it is possible to study the ratio between the second harmonic power 
due to the y vertical and the x horizontal polarization components, that is only 
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Fig. 5. Solid line: illustration of the behavior of F2/F1 as a function of N. Dashed 
line: its asymptotic, Tr/{AN), for N 1. 

a function of the K parameter and is simply given by the ratio R — W2y/W2x' 



R{K) 



2?/ 



w- 



2x 



A\K) 



(58) 



A plot of R{K) is given in Fig. 6. As it is seen the relative magnitude scales 
from 4% in the case -C 1 till about 6% in the limit 3> 1: as one can see, 
the vertical polarization component of the radiation depends quite weakly on 
the K parameter. The knowledge of the polarization contents of the radiation, 
even if relatively small as in this case, can be important from an experimental 
viewpoint. For example, in the VUV wavelength range, the reflection coeffi- 
cients of many materials (e.g. SiC, that is widely used for mirrors) exhibit a 
complicated behavior, and there may be even an order of magnitude difference 
depending on the polarization of the radiation. It should be noted that R{K) 
is independent of the particular model chosen for the beam modulation as it is 
easy to understand inspecting Eq. (32). It is also important to remark the fact 
that the second harmonic radiation from a planar undulator is linearly polar- 
ized, since vertical and horizontal polarization components are characterized 
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Fig. 6. Illustration of the behavior of the ratio between the second harmonic power 
due to the y vertical and the x horizontal polarization components, R{K). 

by the same phase factor. This fact is well-known in Synchrotron Radiation 
theory for a single particle and it is true for any observation angle and any 
harmonic of the radiation from a planar undulator in contrast, for in- 
stance, to the case of bending magnet radiation, when vertical and horizontal 
polarization components exhibit a relative 7r/2 phase shift, indicating circular 
polarization. 

An important comment to what has been done before is needed. We calculated 
the electric field, the angular intensity distribution and the power for the 
second harmonic making a particular assumption about the electron beam 
modulation in Eq. (41). This amounts to consider the modulation wavefront 
orthogonal to the direction of propagation of the beam. The same assumption 
has been implicitly done calculating the first harmonic power (the expression 
in has been used, which does not account for deflection angles). In this 
particular case we have seen that the total power of the second harmonic 
radiation does not depend on the deflection angles rj^ and rjy. In the more 
general situation we find that the second harmonic power can be independent 
of the beam deflection angle (like in the situation treated by us) or can decrease 
due to the presence of extra oscillating factors in / in Eq. (39). On the contrary, 
in jSi] , an increase of the total power is reported due to the presence of deflection 
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angles. 



4 Discussion 



After the presentation of our theory has been given, in this Section we want to 
discuss in a more detailed way the differences between our approach and the 
currently accepted treatment of the problem of second harmonic generation. 

It is worth to begin summarizing the steps which led us to our main re- 
sults. First, we started from the wave equation assuming that the electromag- 
netic sources are given externally by some code calculating the electron beam 
bunching at the second harmonic. Second, after applying the ultrarelativistic 
approximation (7^ S> 1) and the resonance approximation [N^j ^ 1), both 
non-restrictive ones, we solved exactly the wave equation using the Green's 
function method. Third, we calculated the angular distribution of intensity 
assuming a given beam modulation and we derived an expression for the total 
power radiated at the second harmonic by integrating the expression for the 
angular distribution of intensity. Finally we compared the expression for the 
second harmonic power with the analogous expression for the first harmonic. 

In and later on in 0, 0] the ultrarelativistic and the resonance approxima- 
tion were used too, but several steps were performed on the wave equation 
which we find incorrect. First, the gradient term in the source part of the 
wave equation is overlooked. Second, the particles trajectory in the transverse 
X direction is expanded (we will comment on this later on), and because of 
that the Fresnel number N is not identified as the main physical parameter of 
the problem. Third, after these manipulations, the wave equation is not solved 
but, rather, the second harmonic power is estimated in the following way: (a) 
the squared of the (manipulated) source parts of the wave equation for the 
second harmonic is calculated; (b) the squared of the source parts for the first 
or the third 0, 0] harmonic is calculated; (c) the ratio between the square 
of the source parts for the second and either the first f3] or the third jsi, 0| 
is taken. As a result, magnitude of the second harmonic power and polariza- 
tion characteristics are predicted that are in disagreement with what we have 
found. In j^, further notions regarding the case of deflection angle between 
the beam and the undulator axis are introduced. We already expressed our 
critical view on this last conclusion at the end of the previous Section. Let us 
briefly comment on the other points mentioned above by analyzing more in 
detail the approach followed in 0, 0, 0] • 

As has been already said, we find that neglecting the gradient term in the wave 
equation is not correct. As we have seen in Section 2 such term is responsible 
for a contribution to the total intensity for the second harmonic both for the 
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horizontal and for the vertical polarization components and, indeed, it cannot 
be neglected. Doing so would result in any case in an overall incomplete result: 
namely one would obtain only part of the horizontally polarized component 
of the field. 

Going further with the derivation in 0, 0, , the motion of the electrons in 
the x-direction is written as a sum of a fast oscillation due to the undulating 
motion and a slow motion due to the betatron functions. On the ?/-direction 
instead, only the slow motion due to betatron functions is present. The beam 
distribution is then considered as a collection of individual point-particles, i.e. 
a sum of 5-Dirac functions. For the i-th electron one may write 

x'i{z) = Xi{z) + Axi{z) (59) 



and 



y[{z) = yi{z) 



(60) 



where Axi{z) describes the fast oscillation, while Xi{z) and yi{z) describe the 
slow motion. 

All the 5-Dirac in the x coordinate on the right hand side of Maxwell equation 
are subsequently expanded as 

6{xi - Xi{z) - Axi{z)) ^ 6{xi - Xi{z)) - Axi{z)6'{xi - Xi{z)) , (61) 



based on the only assumption that the transverse beam dimension is much 
larger than the wiggling amplitude of the electron motion. It should be noted 
that, based on this assumption, the ratio between the wiggling amphtude 
and the transverse beam dimensions, K/ {'^kwu) <^ 1, is identified as the main 
physical parameter of the theory and is denoted as the coupling strength of the 
second harmonic emission. In contrast with this we have found that the exact 
solution of the wave equation depends on the Fresnel number, but not on the 
coupling strength (see, for instance, Eq. (47)). In this regard it is suggestive 
to write = L^c/{a'^uj) as = [K/{-ik^a)f x tiN^{2 + K'^)/{AK'^). 
As K/{'~fkwa) assumes a fixed value (for instance, much smaller than unity 
according to the assumption above), our Fresnel number can assume any value, 
depending on the number of undulator periods A^^. If, on the one hand, A^ ^ 1 
we have a behavior F2(A^) ~ 1/{AN^) oc A'^ so that F2(A^)/Fi(A^) oc A^^ as 
has been seen in Eq. (57) and therefore the ratio W2/W1 is independent on 
A^ as one can see from Eq. (56). On the other hand, if A^ -C 1 we obtain 
F2{N) ~ \n[l/ {AN'^)] ~ const so that F2(A^)/Fi(A^) is independent on A^^ 
while W2/W1 oc N~^. On the contrary, being based on the coupling strength 
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parameter only, current understandings of the second harmonic mechanism 
predict that W2/W1 is always independent on N^j. 

Finally we find that, from a mathematical viewpoint, the expansion in Eq. 
(61) constitutes an incorrect manipulation of the right hand side of Maxwell 
equation and. To show this, we simply need to consider the mathematical 
structure of the wave equation. For any polarization component we are dealing 
with a differential equation: 

CEx,y{x, y, z) = fx,y{z) J2 ^(^i ~ ^i(^) ~ ^Xi{z))6{yi - yi{z)) , (62) 



with 



.^(V.^.^^l). (63) 

where fx,yiz) is a function containing the appropriate phase factor. This is 
essentially equivalent to our starting equation, Eq. (10), the only differences 
being that the sources are presented in a different way and that, at that stage, 
we had already assumed that the transverse beam dimensions are not smaller 
than the wiggling amplitude of the electron motion. 

The problem with the expansion in Eq. (61) is that Axi = Axi{z) is a function 
of the longitudinal coordinate and that the Green's function for the wave 
equation depends on both longitudinal and transverse coordinates. Let us see 
this point in more detail. If we call with G{zo — z', Xo — x', ijo — y') the Green's 
function of the operator C we have 



00 00 00 



—00 —00 —00 

xU,y{z')5{x' - x{z') - Ax{z'my' - y{z)) , (64) 

that is 

00 

Ex,y{zo-, Xo-, yo) = j dz'G{zo - z\ Xo - x{z') - Ax{z'),yo - viz')) 

—00 

xUyi^') . (65) 

It follows that the expansion of the (5-Dirac in Eq. (61) is mathematically 
equivalent to the expansion of the Green's function G in Axi{z') around Xo — 
Xi{z'); however under the only assumptions Xi{z') 3>| Axi{z') \ and yi{z') 3>| 
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/S.Xi{z') I we cannot expand the Green's function in Axj(z') around Xo — Xi{z'). 
In fact we have that 



G{zo - z\ Xo - Xi{z') - Axi{z'),yo - y{z')) 

G{zo - z,Xo- Xi{z'), yo - miz')) 

dG{zo- z',^,yo-yiiz')) 
-^Xi{z) — , (66) 

because G is simultaneously a function of z not only through Axj(z'), Xi{z') 
and yi{z') but also through Zo — z'. 



5 Conclusions 



In this paper we addressed the mechanism of second harmonic generation in 
Free-Electron Lasers. 

We found that an early treatment of this phenomenon j3] is based on arbitrary 
manipulations of the source term of the wave equation, which describes the 
electrodynamical part of the problem. First, an important part of the source 
term is neglected and, second, an expansion of the particles trajectory in the 
transverse horizontal direction is performed, while we find that there is no 
ground for such a step. Moreover, this leads to the identification of the ratio 
between the amplitude of the electron wiggling motion in the (planar) undu- 
lator and the electron beam transverse size as the main physical parameter, 
while such parameter does not play any role in our theory. The same steps 
were also followed in 0, 0|- After these manipulations, the wave equation is 
not solved but, rather, an estimation of the second harmonic power is given by 
calculating the squared of the manipulated source parts in the wave equation 
for the second and either for the first or the third 0, 0] harmonic and, 
subsequently, taking the ratio between the squared of the second harmonic 
source part and either the squared of the first or the third. Finally in g it is 
introduced the notion of a second harmonic power increasing when a deflection 
angle between the beam trajectory and the undulator z direction is present. 
On the contrary, we find that such power can only decrease or, at most, be 
independent of the deflection angle, depending on how the beam modulation 
is prepared. 

By solving analytically the wave equation with the help of the Green's func- 
tion technique we derived an exact expression for the field of the second har- 
monic emission. We limited ourselves to the steady-state case which is close 
to practice in High-Gain Harmonic Generation (HGHG) schemes but, for the 
rest, we did not make restrictive approximations. This solution of the wave 
equation may therefore be used as a basis for the development of numerical 
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codes dealing with second harmonic emission, which should be using as input 
data the electron beam bunching for the second harmonic as calculated by 
self-consistent FEL codes. 

We found that, in general, the second harmonic field presents both horizon- 
tal and vertical polarization components and that the electric field is linearly 
polarized, while the relative magnitude of the power associated to the verti- 
cal polarization component to that associated to the horizontal polarization 
component is a function of the undulator deflection parameter only. Using 
our result we calculated analytically the directivity diagram and the power 
associated with the second harmonic radiation assuming a particular beam 
modulation case. We expect that these expressions may be useful for cross- 
checking of numerical results. 

In this paper, we presented a theory of the second harmonic generation mech- 
anism in XFELs and pointed out several notions on such mechanism which 
we consider incorrect. In this regard, it should be noted that some of them 
appear to go beyond the subject of harmonic generation itself. In fact we 
have seen that these notions imply an increase of the second harmonic power 
when an angle between the beam direction and the undulator axis is present, 
as if this was a general property depending on Maxwell equations. In other 
words, it looks like the solution of Maxwell equations for any electron beam 
in an undulator would yield a non-trivial dependence on the angle between 
the trajectory and the undulator axis. If so, this conclusion should be valid, 
in particular, for a single particle as well. We find that this is not correct: in 
fact, for a single particle in a undulator the dependence of the electric field on 
the angle between the average direction of the particle and the undulator axis, 
rix^y, is simply related to the chosen reference system, i.e., usually, one with the 
z axis aligned with the undulator axis. A simple rotation of an angle rj^^y to a 
system with the z axis aligned with the electron average velocity woTiId give a 
result independent on such angle. This means that the basic characteristics of 
undulator radiation, in particular the intensity distribution at fixed frequency 
and the spectrum at fixed observation angle O^^y, depend on the combination 
{dx,y ~ Vx,y) only. In other words the presence of an angle rj^^y between the 
electron direction and the undulator axis has the only effect of introducing a 
rotation in the expression of the electric field otherwise leaving unvaried all 
its characteristics, including its resonance frequency. 
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